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. . . We discuss the dynamics of particles in one dimension in potentials that are random both in space 

Cn , and in time. The results are applied to recent optics experiments on Anderson localization, in which 

the transverse spreading of a beam is suppressed by random fluctuations in the refractive index. 
^^ ' If the refractive index fluctuates along the direction of the parajcial propagation of the beam, the 

Cn I localization is destroyed. We analyze this broken localization, in terms of the spectral decomposition 

^—^ . of the potential. When the potential has a discrete spectrum, the spread is controlled by the overlap 

, ^ ' of Chirikov resonances in phase space. As the number of Fourier components is increased, the 

resonances merge into a continuum, which is described by a Fokker-Planck equation. We express 
the diffusion coefficient in terms of the spectral intensity of the potential. For a general class of 
V^ ■ potentials that are commonly used in optics, the solutions of the Fokker-Planck equation exhibit 

anomalous diffusion in phase space, implying that when Anderson localization is broken by temporal 
fluctuations of the potential, the result is transport at a rate similar to a ballistic one or even faster. 
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pH ■ For a class of potentials which arise in some existing realizations of Anderson localization atypical 

behavior is found. 
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I. INTRODUCTION 



It is well known that a particle moving in a spatially disordered time-independent potential can exhibit Anderson 
, localization, and that this is the generic situation in one or two dimensions [l|,|2|. However, it is very difficult to see 

"^ unambiguous evidence for Anderson localization in solids, because of complications due to electron-electron interac- 
Q ! tions. Recently Anderson localization was demonstrated by trapping of light propagating in a paraxial disordered 
Q ■ optical system Q, in the scheme called "transverse localization" f4!|. That experiment had much follow up in optical 
'~~' I setting [5|, |6| , while related experiments were carried out also with cold atoms [7] . In the optics experiments in the 
transverse localization scheme, the role of time is played by distance along the waveguide, and considerable care is 
^ ' required to ensure that the disordered modulation in the refractive index is perfectly uniform along the direction of 
QQ propagation (which we take to be the z-axis). 

CO If the potential fluctuates in time as well as in space, the arguments leading to Anderson localization break down. 

OO Anderson's original paper was entitled 'The absence of diffusion in certain random lattices', hence it might naturally 
^^ be expected that when Anderson localization is destroyed - it will be replaced by diffusive transport. In the optical 
CO setting of transverse localization, the spatially-disordered potential corresponds to small random variations in the 
^^ • refractive index in the xy plane, normal to the propagation direction, whereas the temporal fluctuations correspond 
CN ' to adding longitudinal fluctuations (in the z direction) to the disordered refractive index. That is, the refractive 
index would have random variations both in the xy plane and in the z direction, but at different rates. In this paper 
we consider the effects of such evolving disorder, which as we show, destroys localization, leading to expansion of 
the optical beam, by virtue of evolving disorder. We argue that when Anderson localization is broken, for a generic 
realisation of the optical potential, the result will be super diffusion: expansion of the optical wavepacket at a rate 
similar to a ballistic or even faster. 

Several different mechanisms for the breakdown of Anderson localization due to temporal fluctuations of the po- 
tential have been discussed, all of which predict that localisation is replaced by diffusion. Mott [8] considered the 
effect of phonons at low temperatures, and argued that this gives rise to a diffusive motion of the electrons termed 
'variable-range hopping' conductivity. Mott [9] also considered the effects of an AC electric field, and suggested that 
a resonant interaction dominates the low-frequency response. It has been argued that an alternative limiting proce- 
dure leads to a distinct type of diffusive response at a low- frequency electric field, termed adiabatic transport [lOj . 
Here, we show that, in contrast to previous work predicting diffusion, the response to a time-dependent fluctuating 
spatially disordered potential may result in super diffusion. Specifically, the diffusion constants have a sensitive de- 
pendence upon energy. Also, if the potential is time-dependent, the energy of a particle will not be constant. If the 
diffusion constant is a rapidly increasing function of energy, the response to a time-dependent perturbation can be a 
super-diffusive, ballistic or even super-ballistic. 

Anderson localization is a quantum mechanical effect (at least in cases where the potential is not large enough to 
trap the particles classically). However, as the energy increases the effects of quantum phenomena decrease, and a 



classical analysis becomes appropriate. As we will see in the present work, rapid spreading in configuration space is 
related to an increase of the kinetic energy, justifying the use of the classical (particle) picture. Also, as the energy 
increases, the effects of the potential become a weak perturbation. Accordingly, in order to characterize the asymptotic 
behavior in the long-time limit we consider the behavior of a particle moving classically in a weak disordered potential 
that is also fluctuating in time. 

The random potentials which are prepared in optics [3] and atom optics |7l.lll| experiments are naturally described in 
terms of Fourier series where the expansion coefficients are independent random variables. This motivates representing 
the random potentials using their spectral content. In practice, there are a finite number of Fourier coefficients 
(denoted here by N), but this number may be large. We shall therefore also consider potentials with a continuous 
Fourier transform, approximating the limit as A'^ — >■ cx). 

In section ini we consider the classical dynamics of a particle in one dimensional potentials that are random in both 
space and time, emphasizing the diffusive spread of momentum. In the case where there is a finite number of Fourier 
components, the theory is formulated in terms of Chirikov's resonance overlap criterion |12 - [l4J . In the limit as A^ — >■ oo, 
we show how the Chirikov resonances are related to an expression for a diffusion coefficient D characterizing random 
changes in the momentum p. In accord with earlier investigations [I5l4l8l|. we conclude that, for generic random 
potentials the diffusion coefficient has a universal power-law dependence on momentum p, such that D{p) ^ \p\~^ as 
\p\ — ^ oo. This in turn implies that asymptotically in time the average momentum satisfies (p^) ^ t^'^ . The average 
displacement satisfies (a;^) ~ t^^/^ for one-dimensional systems [1&-18] (faster than ballistic) and (x^) ^ t^ (ballistic 
transport on average) for systems with dimension higher than one [16]. 

Our discussion in section |ll] differs from earlier works analyzing anomalous diffusion in random potentials |l5l4l8( , 
in that we express the diffusion coefficient 13 (p) in terms of the spectral intensity of the potential, as well as making 
the connection with the Chirikov resonances explicit. This approach highlights some of the subtleties which can arise 
when considering real experiments. In Section Hill the theory is used to derive the transport properties for a potential 
that naturally appears in optical experiments on Anderson localisation, such as p|. We find a family of potentials 
for which the diffusion coefficient of the momentum vanishes for high momentum. This suggests that the diffusive 
spreading in momentum saturates asymptotically, and therefore does not exhibit the universal behavior described in 
previous studies (l5l - [l8| . In the optical experiments such as [3)SQ both types of potentials could be readily realized, 
which allows the exploration of both transport regimes. The formulae presented here give a quantitative prediction 
of the anomalous spread of the beam. The results are summarized in Section lIVl 



II. A PARTICLE IN A QUASI-PERIODIC POTENTIAL 

Some of the experiments which demonstrate an optical realization of Anderson localization (such as that described 
by [3]) involve an induction technique, where a change in the refractive index of a dielectric is induced by an interference 
pattern generated by external waves (used strictly to induce the potentials) [19, 2Q|]. Because the optical field defining 
the disordere is produced by interference of plane waves, the 'potential' is most naturally described in terms of its 
Fourier components. For this reason we need to analyze motion in a quasi-periodic potential. 

Consider the motion of a particle of unit mass described by the Hamiltonian, 

il=t^^V{x,t\ (1) 



where Vix^ t) is a one-dimensional quasi-periodic potential of the form, 

N 

'N 



1 ^ 
V{x,t) ^ —j= ^ A„iex]i[i{k„iX-u„rt)] (2) 



= -N 



where A-m = ^mi ^^ that the potential is real. Here Am are independent (for to > 0), identically distributed complex 
random variables. The expectation values of these variables satisfy (for n,m > 0) 



\Am) — \AmAn) — U 



{AmAD = a'^Smn. (3) 



An example of such a variable is Am = \Am\ exp (i0m), where \Am\ and </>,„ are independent real random variables, 
with (j)m uniformly distributed in the interval [— tt, tt]. The random variables km and ujm are distributed with the 
probability density P(fc,w), which may be either a continuous spectrum, a sum of delta functions, or a distribution 



concentrated on a line in k-u space. Note that, 

N 



\Vix,t)\< V |A„|/ViV<2\/7Vmax|An|, (4) 



m=-N 

however the variance of |F(a;,i)| is equal to 

N N 
N 



\V{x,t)\'^j = ^^ X! ^{ArnA*^)exp[i{{k.m-kn)x~{uJrn-(^n)t)] 
m.= — N n=l 
N N 
— ^ ^ (AmAn) exp [i ((fc„i + kn)x - {UJ„1 + UJn)t)] . (5) 

Using the assumptions of ([3]), we have 



N 

m— — N n— 1 



\V{x,t)\')^2a\ (6) 

For finite N the potential is a quasi-periodic function of x and t, and in the limit of A'^ — > oo, for fixed x and t it is a 
Gaussian random variable with the variance, a^ . 

The motion of a particle in a potential given by ([2|), for sufficiently small \Am\, (the exact requirement will be 
specified below) was analyzed by Chirikov [13]. It was predicted that the phase-space is built up of chains of non- 
overlapping resonances, which are given by the condition, 

— (k„iX - UJmt) = k,-nP - UJrn = 0, (7) 

at 
which is just the stationary phase requirement. This reduces to the condition, 

P™ = ^. (8) 

Assuming that the resonances are isolated, starting a particle with an initial momentum near a resonance, will 
produce a bounded pendulum-like motion near that resonance. This can be seen by neglecting all non-resonant terms 
in the potential, and making a Galilean transformation to the frame of reference of the specific resonance. The new 
Hamiltonian in this frame of reference is just the time-independent Hamiltonian of a pendulum, 

H'^^^+'^dcosk^x. (9) 



We can estimate the width of the resonances. A™ (that is, the range of momentum for which phase points lie on 
oscillating trajectories). From energy conservation, A^/2 = A\Am\l^fN, 



\J^a„,\/Vn. (10) 



We now order the resonances, such that, p^i^ < p'2^ < ■ ■■ < jp^y, and define the distance between the adjacent 
resonances by, Sm = Pm^ — Pm-i- The Chirikov criterion |12 - ll4l | for trajectories to remain localised close to their 
initial momentum is that the resonances do not overlap, that is 

(A™-fA„._i) < S,n Vm, (11) 

Under these conditions, the momentum of a particle will not change appreciably over time. One should note that 
condition (TTl]) is approximate. To obtain better estimates higher order resonances should be considered (l4| . 

When the amplitudes of the potential, \A„i\, are not sufficiently small, namely. Am ~ Sm, some of the resonance 
chains will overlap. It is established [12, 13] that in places where resonances overlap stochastic regions will form, which 
will result in a random walk between resonances and therefore a diffusion in the momentum. In order to observe 
diffusion the number of resonances has to be large, since for the diffusion approximation to be valid a large number 
of jumps between the resonances has to occur. We will now obtain the diffusion coefficient, adapting a technique 
developed in [l5l - [lg ] to the case where the potential is described by the statistics of its Fourier components. In the 
limit as N ^ 00 the Chirikov resonances become dense in momentum, which appears at first sight to complicate 



the problem. However, in this hmit the quasi-periodic potential is replaced by a random potential, and the change 
of the momentum in a time interval which is longer than the correlation time of this potential can be regarded as 
a stochastic variable. In this limit, the small changes in momentum which occur over a timescale which is large 
compared to the correlation time of the potential can be treated using a Markovian approximation, which validates 
the use of a Fokker-Planck approach. 

We will now proceed in line with [13], writing an expression for the small change in momentum occurring over a 
time St which is large compared to the correlation time of the potential 

pSt 

5p= dt F{x{t),t) (12) 



where x{t) is the trajectory of the particle and F{x,t) = ^ix,t) is the force. Defining the force- force correlation 
function, 

Cixuh; X2M) = {F{xi,ti)F{x2,h)} , (13) 

and assuming that it is stationary, we can express the variance of the fluctuation of the momentum in the form 

{Sp^)= dh dt2C{x{h)~x{t2),h-t2), (14) 

Jo Jo 

and we neglect all the cross-correlations {SpiSpj)^ where the indexes i and j correspond to two different intervals. 
For this assumption to be true the correlation function C should decay sufficiently fast, such that for t > St it is 
negligible. Furthermore, we will expand 

x{ti) - X{t2) = p{ti - t2) + 0{St^), (15) 

which assumes that the force and its time variations are weak enough. Under these assumptions, we can obtain, 

(V) = ^Dip)St, (16) 

where the diffusion coefficient is given by, 

D{p) = ^J^ C{pT,T)dT. (17) 

This expression was first obtained by Sturrock [21| . As a consequence of the increments of momentum being both 
small and Markovian, the probability density for the momentum, p{p,t), satisfies a Fokker-Planck equation. Care 
must be taken over the order in which derivatives are taken. In [18i| it is shown that the correct Fokker-Planck equation 
is 



Using the potential ^ we obtain the correlation function, 

N N 

N 



N N 

C{xi,ti]X2,t2) ^ — ^ ^ (fc„fc„A„exp[i(fc„a;i -a;„iii)] X y4* exp[-i(A:„a;2 - w„t2)]) • (19) 



m=~N n=-N 



We first perform an average on the A„i variables. Using the assumptions (j3]) results in a translationally invariant 
correlation function both in space and time, 

2 ^ 
C{xi-X2,ti-t2) = -^ X! (^mexp[i(fc„i(a;i -2:2) -a;„i(ii - ^2))] +C.C.) 

m=-N 

= a"^ j dk I duj k'^P{k, oj) (exp[i(fc {xi - X2) - oj{ti - t2))] + c.c.) , (20) 

where P{k,uj) is the probability density of uj and fc, which will dubbed in what follows the spectral content of the 
potential, introduced along with equation ([3]) . Note that the correlation function of the force is a Fourier transform 
of k^ P{k,uj). To obtain an integrable correlation function, we therefore have the following requirement on P{k,uj): 

f dk f dujk^P{k,uj) = M <oo, (21) 



which means that P{k,uj) has a finite support or decays faster than k ^ and ui ^. Note, that distribution P{k,uj) can 
contain also atom contributions (delta functions). Equation ([20)1 leads to, 

D{p) = ^ dk duj fc2 P(fc, w) / dr (exp[i(fcp - lu)t] + c.c.) 

= 2tt(j^ dk dujk^ P(k,uj)S{uj-kp). (22) 

This expression is closely connected to the resonance probability density, which can be defined as 

P(p-^) = I dk I dujP{k,uj)5{p'^'-'^) (23) 

= f dk f duj\k\P{k,uj)S{uj-kp"'''). 



These two expressions quantify the connection between Chirikov resonances and phase space diffusion. 

Equation (j22p can be also interpreted as an integration of the function fc^ P{k, w) over a line with a slope of p, which 
is just a Radon transform. We can obtain the asymptotic behavior of (j22p for large p following a similar procedure 
done in (l7l.ll8l|. by rescaling the variables, k' = kp, 



9 2 r<yo y 

D{p)^^ dk'k'^P{-X)- 

Therefore if P(0, k') 7^ in the limit of large p we have. 



(24) 



where, 



Dip) ^ ^, (25) 

/OO 
dk'k'^P{0,k'). (26) 

-00 

This scaling of the diffusion coefEcient provides an anomalous diffusion in momentum such that (p^) ~ t^^^ and 
(x^) ~ i^^/^ (l5l.ll6fl. The precise prefactors can be extracted from results in [l7i| (equation (30)) 

^^.^^^'^'Mn/mmm^)^v.^i^ (27) 

The prefactor in the relation 

{x^)^Cj^. (28) 

is given in flq in terms of a one-dimensional integral (equations (150)-(152)). 

In our discussion of optical realisations of Anderson localisation we will be led to consider potentials for which there 
is some pmax, such that 

P{k,Lu) = : w>p„axfc, (29) 

namely, the function P{k,ui) has a non-vanishing support only outside the wedge with an intercept of pmax, than it 
fohows from ([22|) that. 

Dip) = h""^^ -^0°° '^^ ^^^ ^^'^^^ '^' - ^"^'^ (30) 

[0 \p\ >_Pmax- 

This suggests that under the condition ([29)) on the potential ([2]) there is a saturation in the growth of the kinetic 
energy. Therefore a particle started inside the part of space with non-vanishing resonance density will not diffuse to 
regions of zero resonance density. 

To summarise: we have derived the explicit dependence of the diffusion coefficient on the spectral content of the 
potential, P(k,uj), which is very useful since in many cases in optics and in atom optics it can be experimentally 
controlled. An example of this type will be discussed in the next section. 



III. APPLICATIONS TO OPTICS 

In this section we will apply the general scheme for the calculation of the diffusion coefficient to a specific realization 
of a disordered potential. In recent experiments examining Anderson localization of light, the potential was realized 
by a superposition of plane waves, which is of a structure similar to ([5]) ^i^]. In those experiments light propagates 
paraxially in a disordered potential: the signature of localization is that the width of the propagating beam of light 
from a coherent and a monochromatic source remains bounded as it propagates. The disordered potential is produced 
by utilizing the photosensitivity of the medium. A powerful polarised writing beam induces a change in the refractive 
index. An, of the medium. The polarization of this beam is selected in such a manner that the beam does not 
experience the change in the refractive index that it induces. The localization experiments are carried out with 
another beam (probe) with a different polarization, such that it experiences the written change in the refractive index 

where Eq^t) is the magnitude of the electric field of the writing beam at position r, and _B is a coefficient proportional 
to the nonlinear susceptibility of the medium [3|, |6| . Iq is some constant background intensity. In this work we will 
consider that \Eo{r)\^ /Iq is small (as is in fs'l) so that we can expand, 

A„„ = i,(l-M£)!!), ,32, 

furthermore we will assume that the refractive index depends upon just two coordinates x and z, where z measures 
distance along the axis of the test beam. We extend the analysis of last section for a particular potential realization, 
which is used in many experiments in optics, showing that when Anderson localization is broken, and the spread of 
the beam obeys an anomalous diffusion law. We provide numerical results to support our conclusions. 

The propagation of a monochromatic light beam in a medium with a non-uniform refractive index is described by 
the Helmholtz equation for the electric field E, 

V^E + e^^E^O, (33) 

where k = l-Kn^jX, A is the wavelength of light in vacuum, no is the bulk refractive index, and n(v) = tiq + An(r). 
Setting, 

E{x,y,z)^i,{x,y,z)e''\ (34) 

for A?i(r)/no ^ I and when 



d^ip 



dz^ 



<2fc 



dz 



the paraxial approximation is invoked. This approximation yields a Schrodinger like equation for the slowly varying 
amplitude, ip {x, y, z). 



-d^ip = -— ^ {dxx + dyy) il) 

k 2k^ no 



dz^ = -7n:2^^^^^^Vv'>'^~ ~~'^' (^^) 



where we have chosen z to be the propagation axis. In this work we will consider a propagation along a one-dimensional 
potential, such that motion is confined to a plane, with the y degree of freedom frozen. 

The fluctuations of the refractive index are achieved by utilizing the sensitivity of the medium to light at some 
frequency and polarization, which allows us to transform a pre-designed interference pattern to a variation in the 
refractive index An{x,z). The interference of N plane waves induces a fluctuation A7i(r) which is equivalent to a 
potential of the form of ^ . The total electric field of the writing beam is given by, 

1 ^ 
Eo{x,z) = -r= ^ Amexp[i{kx^jnX + k^^rn)z] (36) 

V ^ m=-N 

with A-jri = ^mi where kx^m and kz^m are the x and z components of the wave- number of a plane wave labelled by 
an index m, for which the magnitude of the wavenumber is fco- The normalization is chosen such that the total power 



does not change as a function of N, 

N 



m=l 



(37) 



where V is the volume of the system, does not depend upon N. The resulting change of the refractive index is equal 
to 

B ^ 

An{x,z) ^ ^\Eo{x,z)\'^ "^^ TTT T^ Aji^^exp [i ((fc^,,™ - fca;,„)a; + (fc^.™ - fc^.n)^)] . (38) 

l\Io '-^ 

The experiments are typically performed with the wavevectors of the driving electric field Eo{r) close to the z-axis. 
This implies that we can use the paraxial approximation for the writing field i?o(r), as well as for the weak probe 
field. This justifies the following additional paraxial approximation: 



2ko ' 
to simplify the notation we will set 



kz,m — \/ ^Q fca;,m ~ '^O r,i ' (^9) 



A-2 

k = k I.) - ^'"' (Af)] 

2ko 



B ^ 



Then ((38)) simplifies to 

An(a;, z) = -— - y^ A„A* exp [i ((fc™ - /c„)a; - (w™ - a;„)z)] . (41) 

JVio „ 

To simplify comparison with the previous section, we will work in units where fco = 1 a-nd /q = 1, and will designate 
the paraxial axis by t (instead of z). Using these conventions ([35)) simplifies to, 

idtE^-]^d^^E + V{x,t)E, (42) 

which has the form of the Schrodinger equation in one dimension. In our model we take the potential to be 

AT 

TV 



1 ^ 
^(2^7^) = T^ X! ^m^^exp[i((A;„j - A;„)a;- (w„-a;„)i)] , (43) 



where we have absorbed all the constants, including the minus sign inside the A^, and where we have confined to 
positive indices because this expression is automatically real. Hence, (|42p describes a propagation of a paraxial light 
beam inside a medium with spatially- varying refractive index. 

Having defined the 'potential' function for the paraxial equation by (|^5)) . we now consider its ray dynamics. The 
classical (zero wavelength) system corresponding to (|42)) is a 'particle' moving in the potential given by (|43)) . The 
equations of motion of such 'particle' are given by 

dp dV , , da; , , , , 

where p is the conjugate momentum to x. In what follows we will examine the motion of such a 'particle'. We will 
now proceed with a similar analysis to the one done in the previous section. Comparing ([2]) and ()43p , we see that the 
resonances of the system are given by 

— [(km - kn)x - (ujm " W„)t] = (fc,^ - fc„)p - {uJ„i ~ UJ^) = 0. (45) 



res ^ ^rn ^n ^ -rn . ■"■n , .„. 



Or using the definition (|40 
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Figure 1: (color online) A phase-space portrait for three resonances. Different colors (shades) describe different trajectories. 
The horizontal black solid and dashed lines, designate the edges of the resonance chains. The initial conditions were uniformly 
distributed on the p axis {x{t — 0) ~ 0). The system was integrated up-to t = 10* (dimensionless variables). At the left figure 
three non-overlapping resonances are presented, while in the right figure resonances (1) and (2) overlap. 



The number of resonances is N{N — l)/2, and their width is, 



A„„ = y/8\A^An\/N. 



(47) 



The overlap condition of the resonances is the same as in the last section. 

The change of 'momentum' (that is, angle) for a ray which propagates in between the resonances or outside of the 
region of resonant chains will average out to zero. This is illustrated in figure [1] where one can see a transformed 
phase-space of such a system with three non-overlapping resonances. The phase-space was sheared by defining 
x'{t) = x{t) — p{t — 0)t. In this manner, trajectories for which the momentum does not change considerably during 
the motion of the particle, are slowly evolving with respect to the transformed x axis, whereas trajectories which 
exhibit a large variation in the momentum show a large excursion in x' . From figure [1] we see that trajectories started 
within a resonance chain show a large excursion in x' . Nevertheless, all of them are bounded by the boundaries of the 
resonance chain. Trajectories started in between the resonance chains or outside of them have a rather small variation 
in the momentum, which decays with increasing the absolute initial momentum. 

In figure [1] (right) we illustrate an overlap between two of three resonances of a system plotted in Fig. [T] (left). A 
transition between the overlapping resonance chains is clearly seen. 

As the number of resonances is increased, we can reach a regime where a large number of these resonances overlap. 
In this regime we anticipate that the propagation angle of the ray, p, could exhibit anomalous diffusion. For the 
regime when all the resonances overlap, we will compute the diffusion coefficient. The correlation function for the 
force is 

N N 

C{xi,tuX2,t2) ^ ^ Y. Y. (^™^«^*^j) 



n,m—l ^,J — 1 



X {{km - kn){ki - kj)exp [i {{km - kn)xi - {uJm - UJn)tl ~ [h - kj)x2 + (a;,j - UJj)t2)]) ■ (48) 



Similarly to the previous section we assume 



\A-mAj^Aj^Aj) — a {OnmOij + OmiOnj)- 



(49) 



This reduces equation (H5)) to 



C{Xi- X2,ti -t2) = -r^ X! {(km- kn)^exp[i{{k„i- kn){xi -X2) - {uJm-i^n){tl-t2))]) 

n,'m—l 

= a^ fdki j dk2P{ki)P{k2){ki 



k2? 



X exp 



i ( (fcl - k2){xi - X2) -\{kl~ kl){ti - t2) 



(50) 



where P{k) is a density of resonances in k space. We used the fact that k and w are related by a dispersion relation 
(|40| . Using the expression for the diffusion coefhcient (ITTI) , we obtain 

Dip) = ^J C{pT, r)dT = Tra^ /" dfci f dfca P(fci)P(fc2)(fci - fc2)'(5 ('(fci - fca)^ - ^(fc? - ^2) 



2' 



,5(fci-(2p-A:2)) 



^ 7:a^ Jdk,Jdk2P{k,)P{k2){ki-k2)^- 

= iira^ f dk2 P{2p - k2)P{k2)\k2 - p\. (51) 

Note that if P{k) has a finite support, than D{j)) will also have a finite support. For example for P{k) = 9{k]i — 
\k\)/{2kfj), where 0{x) is a step function, using ([5T|) we compute the diffusion coefficient 

Notice that this expression has finite support, so that the asymptotic relation which leads usually to anomalous 
diffusion, ([25]) , is not satisfied. In this case a particular form for the potential, which is used in optical realisations of 
localisation such as |[3|], leads to sub-diffusive growth of (p^). For the case of a large amplitude of the potential, where 
all resonances overlap, the variation of the momentum of the trajectories is bounded, in accord with the prediction 
that D{p) has finite support. 

In figure[5]we compare a direct numerical evaluation of (p^(i)), obtained by averaging solutions of pi)) over different 
initial conditions, with an estimate obtained from a numerical solution of the Fokker-Planck equation 

which uses ([5^ and k^ — 0.1. The classical trajectories all had initial conditions in the resonance chain. A good 
correspondence is found without any fitting parameter. 

IV. SUMMARY AND DISCUSSION 

We have argued that Anderson localization may not be the absence of diffusion in disordered potentials, but rather 
the absence of anomalous diffusion. More precisely, using the standard argument that at high energies quantum 
interference effects are negligible, we used effective particles to describe the dynamics, and argued that when Anderson 
localization is destroyed (due to temporal fiuctuations of the potential) the long-time dynamics is determined by a 
semi-classical approximation. For generic potentials, this semi-classical dynamics exhibits anomalous diffusion in the 
long-time limit. 

Anomalous diffusion in a disordered classical model has previously been studied in several works |l5l4l8| . This work 
has considered the issues which arise when considering whether this type of anomalous diffusion is also relevant to 
the breakdown of Anderson localization in optical systems. 

The context in which Anderson localization is most readily accessible to experiment is in propagation in a disordered 
optical potential (refractive index) [S] induced by an optical interference pattern. For this reason, we concentrate upon 
the case where the potential is quasi-periodic, resulting from the addition of N waves. We discussed the influence 
of the potential on paraxial propagation of a coherent beam, showing how a semi-classical analysis is related to the 
Chirikov resonance overlap criterion |12l . Il3l | of Hamiltonian dynamics. 
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Figure 2: (color online) (p^) as a function of t. Dashed gray lines are numerical solution of (|44[) for 1000 initial conditions taken 
from a Gaussian distribution. Different shades and styles of lines designate different A'' (see legend). The solid think black lines 
are the numerical solutions of (|53|) with a diffusion constant given by, (|52p . The parameters are, ku — Q.\ and a = 10""' for 
the left figure and a = 10~* for the right figure. 



Gradually increasing the amplitude it is demonstrated how one transforms from the regime of isolated resonances 
where no diffusion takes place (Fig. [TJ left) to a situation where few resonances overlap (Fig. [U right) and spreading 
that involves them is found. When the amplitude is even further increased, such that all the neighboring in momentum 
resonances overlap and diffusion in momentum is found (Fig. [5]). 

In the limit as the number N of Fourier components increases, the potential is described in terms of a spectral 
intensity function P(fc,w), and the effect of the potential can be modeled by a stochastic equation, describing fluc- 
tuations of the momentum p with a diffusion coefficient D{p). We showed how D{p) can be related to the spectral 
content P{k,uj) by ([^^ . and how this relation can be interpreted in terms of the Chirikov resonance condition. We 
also showed that earlier results on anomalous diffusion, ([?7)) and (|28l) are recovered. 

We also considered in some depth the type of potential, such as (|^5| . which arises in optical realizations of Anderson 
localization such as [3|. We showed that the potential which is used in these experiments is non-generic, and leads to 
(p^(t)) being practically bounded, rather than exhibiting anomalous diffusion. This implies that experiments to test 
the prediction that breaking localisation leads to anomalous diffusion will have to be carefully designed, such that 
the span of the momentum spectrum of the spatial disorder exceeds the plane- wave spectrum of the initially bounded 
probe beam, in order to observe anomalous transport. 
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